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We derive an effective classical theory for real-time SU(iV) gauge theories at high temperature. 
By separating off and integrating out quantum fluctuations we obtain a 3D classical path integral 
over the initial fields and conjugate momenta. The leading hard mode contribution is incorporated 
in the equations of motion for the classical fields. This yields the gauge invariant hard thermal loop 
(HTL) effective equation of motion. No gauge-variant terms are generated as in treatments with 
an intermediate momentum cut-off. Quantum corrections to classical correlation functions can be 
calculated perturbatively. The 4D renormalizability of the theory ensures that the 4D counterterms 
are sufficient to render the theory finite. The HTL contributions of the quantum fluctuations provide 
the counterterms for the linear divergences in the classical theory. 



I. INTRODUCTION 



To understand the processes that have played a role in the early universe, it is important to study the behavior 
of thermal field theory at high temperatures. Some of the processes with cosmological implications, like sphaleron 
transitions and the dynamics of weak first-order transitions, are sensitive to the soft modes (k ~ g 2 T, with g the 
gauge coupling) of the magnetic sector in hot gauge theories [j^] . This low- momentum behavior is non-perturbative, 
due to the IR problems of massless field theories in 3D Jj]|| . 

For static magnetic fluctuations the dominant IR behavior can be isolated from the full quantum theory by the 
technique of dimensional reduction ^J^] , which allows one to study the non-perturbative behaviour by standard lattice 
simulations or other imaginary-time methods HQ . In a number of papers it has already been argued that real-time 
thermal field theory in the high-T limit reduces to a classical statistical theory [0 Such a theory consists of equations 
of motion that govern the time-evolution of the classical fields and an averaging over initial fields with a thermal 
weight. The goal of this paper is to obtain an effective classical theory by integrating out quantum fluctuations using 
the method similar to that of dimensional reduction || . 

A major problem for gauge theories is that the perturbative fluctuations (hard modes) do not cleanly decouple from 
the soft modes one is interested in. Even in lowest order one will have to take the hard-thermal-loop (HTL) effects 
into account . To handle this in a gauge invariant manner we will propose in the present paper a way of splitting 
the field in perturbative and non-perturbative modes, other than the standard method of introducing an intermediate 
scale to separate the soft from the hard modes jTo|— p~3[| . In this way the influence action obtained by integrating out 
the hard modes is BRS invariant, in the HTL approximation this implies gauge invariance [fl4|| . 

the theory that we obtain after the classical approximation is very similar to the effective theory proposed by Iancu 
|l3"|| . It consists of the HTL equation of motion and a thermal average over initial conditions. The approach that is 
taken here shows that the HTL's should be calculated with a subtracted propagator. The subtraction follows directly 
from separating off the path integration over initial fields from the quantum path integral. This subtraction is essential 
to justify the classical approximation, and it introduces automatically the counterterms for the linear divergences of 
the classical theory that otherwise have to be introduced by a rather ad hoc matching procedure |D| . 

The outline of this paper is as follows. In section [n] we separate off the quantum fluctuations and define an 
influence action containing the quantum corrections to the classical action. We give the Feynman rules to calculate 
this action and use it to define the thermal weight and effective equations of motion of the effective classical theory. 
In section III we calculate the leading order quantum corrections to the thermal weight. We consider especially the 
terms {£ l f~£ l V l A°, (V l A°f and (A ) 2 , since their coefficients may be different in the Hamiltonian and Lagrangian 
formulations of static dimensional reduction. The leading order contribution is the Debye mass term, the only static 
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HTL. We then argue that the HTL contributions should be included in the effective equations of motion. The HTL's 
obtained through integrating out the quantum fluctuations are linearly divergent. These are, in fact, the counterterms 
for the linear divergences encountered in the effective classical theory. In section IV we present the effective classical 
theory and discuss the classical approximation. Section |v| contains the conclusion. 



II. SEPARATING OFF THE QUANTUM FLUCTUATIONS 

A. Background fields and zero modes 

In the real-time formulation of thermal field theory, the thermal state is taken into account by extending the time 
evolution along a contour C in the complex t-plane, as shown in fig. 1. The contour consists of three parts: a forward 
time branch C\ on the real axis starting at an initial time t m and ending at a final time U, a backward branch C2 
from t{ to t m and a third branch C3 down the Euclidean path t = t- m — it from r = to r = (3. On this contour we 
consider a SU(iV) gauge theory for the field = (A , —A) with generating functional [^5|Jlq|, 



Z[j] = J DADEDcDc <5(V • A) exp (iS[A, E] + iS gh [c, c, A] + ij ■ A^j , (1) 

where E is the momentum conjugate to A. The ^-function enforces the Coulomb gauge V • A = 0. We use the 
inner-product notation j ■ A = —2trf c d 4: xj fJ/ (x)A ,i (x); the notational conventions are those of reference pjjfl . The 
action of the theory is 



S[A, E] = 2tr / <Px 



4 ■ 

c 



E.d t A+-(E 2 ! If) \"n 



(2) 



Here B is the magnetic field and A is the Lagrange multiplier associated with the constraint R(A, E) := V • E + 
g[A, E] = 0; it is real on the real-time contour and complex on the imaginary time- branch A°(t- ln — it, x) = iA (r, x). 
The time-derivative of the field A in (||) is taken along the contour. The integration over the ghost action 

Sgh = -2tr / d 4 x (cV 2 c - gc [A, Vc]j (3) 

generates the Faddeev-Popov determinant. The fields A and the ghosts are periodic with respect to begin- and 
endpoint of the contour, but the momenta E and the Lagrange multiplier A have no boundary conditions. 

We introduce the classical background fields on the real-time contour Cyi — G\ U C2 by shifting the fields A —> 
A c i + A, E — > E c \ + E, and similarly for the ghost fields. Performing this shift in the action (|2|), we obtain 

S[A cl + A, E c i +E] = S[A, E] cl + 6 A S[A, E] cl ■ A + S[A ch E cl ; A, E] , (4) 

where we have defined S[A, E] c \ = S[A, E]\ A=A j E=E . The last term in (||) is the part of the shifted action containing 
terms quadratic-and-higher in the quantum fluctuations. The term linear in E has been made to vanish by imposing 
on C12 the condition 5eS[A, E] c \ = which defines the classical electric field in term of the potentials 

E l cl = d l A°rt d°A l cl g[A l cU A cl } =: F<° . (5) 



on C\2- In section [V, we will require the background field A c \ to satisfy the hard thermal loop (HTL) equations of 
motion. 

The procedure for treating the ghost action <j3j> is the same as for the gauge field action (^). Linear terms are 
eliminated by imposing S s S g h[c, c, A] c \ — 0, S c Sgh[c,c, A] c i = 0. Since the ghost fields are not dynamical in the 
Coulomb gauge, this implies c c \ — 0, c c \ = 0. 

On the Euclidean part of the contour we follow the prescription of static dimensional reduction and split the fields 
into zero and non-static modes. The zero modes are constant in time and defined by a projection of the fields on C3; 
for example 

A(x) = VA(t, x)=iT dtA(t, x) . (6) 
Jc 3 
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The other zero modes are defined similarly and denoted as £,C,C. The static and the non-static modes are separated 
in the integration measure by formally writing DA — > DA [DAS (VA — A)] , and similarly for the other fields. Since 
in the integration over the quantum fields A the zero mode A is treated as a constant background field, we may shift 
A — > A + A like on the contour Cn. The shifted action has the same form as in (Q) if we identify A c \(x) — A(x), 
E c \{x) — £(x), C c \(x) — C(x) and C c \(x) — C(x) on C3, except that there is no linear term. This is a well known 
feature of static dimensional reduction. In the present case one may verify that the linear terms vanish on C3, because 
they contain the projections VA, etc. 



B. Influence action 

Performing the manipulations described in the previous section, we separate the generating functional ([l]) into a 
path integral over static fields 

Z[j] = J DAB£DCDC5(V ■ A) 

exp (iS[A c i, E c i] + iS g h[cch c c i, A c i] + iW[A c i, E ch c c i, c c i; J] + ij ■ A c ij , (7) 
and one over the quantum fields which defines the influence action 

exp iW[A ch E ch Cci, c cl ; J]= J DA^DcDcS (VA) 8 (Vc) 5 (Vc) <5(V ■ A) 

exp (iS[A ch E cl ; A] + iS gh [c c i,c ch A c i;c,c, A] Aj . (8) 

Here we used that the gauge constraint may be imposed separately on the zero and non-static modes. The source J 
is the combination of the external source and the linear term in (Q) on Cyi 

J[i = 3^ + ^A^Sci . (9) 

The ghost action 5* g h[c c i, c c \, A c \; c, c, A] contains the terms quadratic- and- higher in the quantum fluctuations arising 
from the shifted ghost action. The Gaussian integration over the momenta E has already been performed, which 
reduces the action in the exponent of (d) to the expression 

S[A ch E cl ;A]=-2tvJ^d 4 x (-\f^F» v + l -g[A k , A 1 ] F$ - gA° [A 1 , E l cl ] - Q^j , (10) 

with the classical field strength F c f = d k A l cl — d l A k } — g[A k l} A l cl ], k,l — 1, 2, 3. The quantum field strength tensor 

= d^A v - v M l - g[A% A v \ - g[A», A v cl ] - g[A», A"] , (11) 
contains linear contributions from the background field. The last term in ( |l"o| ) is given by 

Q = l( VF oif } (12) 

arises from the integration of the momenta and lives only on C3. 

We note that on the real-time contour C12 the action ([To]) is just the standard quantum action for SU(iV) shifted 
by a classical background: 

S[A ch E cl ; A] + J ■ A - S[A cl + A] - S[A cl ] +j-A. (13) 

To derive ([l3]) we inserted the equations of motion (||) for E c \ into the action (Q), which gives the Lagrangian version 
S[A] := S[A, E — F 10 ] for the SU(A r ) action. On C3 the electric field E c \ — £ is an independent integration variable 
and we cannot simplify to a form like (p_3J). 
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C. Feynman rules 



The effective action can be calculated perturbatively with the Feynman rules determined by the action (|l^). By 
inspection of ( [j"3| ) we see that on C12 the interaction vertices are the usual vertices with two, three or four quantum 
fields on the external lines. Hence, the Feynman rules for a diagrammatic evaluation of vertex functions on C12 consist 
of the usual interactions. But the propagator of the quantum fluctuations is not the usual propagator of thermal field 
theory, as we will derive below. 

On C3 there is an interaction containing the static electric field, this interaction gives corrections, for instance, to 
the £ 2 term in the classical action. There are also interactions contained in Q (|lj) these do not contribute at one-loop 
order and are not considered in this paper. However the quadratic part of Q is however essential in the computation 
of the propagator of the quantum fluctuations. 

The influence action (ph for free fields only, depends on the background fields through the source (g) and will be 
denoted by Wq [J] ■ This quantity determines the propagator A 1 * 11 associated with the quantum fluctuations according 
to 

ex P iW [J]=e X p--j dtdt' J j—^ J„(t, -k)Ag(t - t', k) J v (t', k) , (14) 

where the tilde indicates the spatial Fourier transform. To evaluate the left-hand-side we insert an integral representa- 
tion S (PA) = J T>Xfj, exp (— x ■ A) for the delta-function in (||). The auxiliary field x is zer0 everywhere except on C3 
where it has a spatial dependence: x^(^ ®) = Xni^) j t € C*3- The quadratic part of the term Q, equation (12), reduces 
to Qo — {VdtA) 2 , since VA vanishes on account of the delta- function constraints in (||). These manipulations allow 
us to rewrite the free effective action in the form 



expiW [J] = exp (^VdtSj ■ Vd t S.^j J D cXai exp -i( J + i x ) ■ D 



c-(J + ix), (15) 

with the standard thermal propagator of gauge fields in the Coulomb gauge JTtJ . The contribution of the Q term 
has been taken out of the functional integration by replacing the spatial gauge fields by the appropriate functional 
derivatives. 

The propagator A^, has the same tensor structure as the standard thermal propagator in the Coulomb gauge: 

A^(t - t', k) = T" v (k)A c (t - t\ k) + g^g u0 A c (t - t', k) . (16) 
The matrix T is the transversal projection operator with Tq v = T^q = and 

T a = % - . ( 17 ) 

with k — \k\. The evaluation of the temporal part of the propagator ( |l6| ) is simple, since the functional derivatives in 
( |i"5| ) do not contribute and the thermal propagator &q is local in time: 

A°S{t - t', k) = D° c °{t - t', k) - S^{t - t', k) , (18) 

with subtraction 

The calculation of the scalar propagator Ac for the spatial gauge fields is more involved; we refer to || for details. 
The result is 

A c (t - t', k) - D c (t - t', k) - S c {t - t', k), (20) 
with L>c the standard thermal propagator [0. The subtraction is given by || 

|j- cos fc(i-i') t,t'eC u , 

^'-''^<|3-y < 2i) 

TZ t,t G C3. 
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On the real-time part of the contour we recognize the subtraction as the classical propagator. On the vertical part 
of the contour it is the zero-mode propagator, which implies that static fluctuations on C3 do not contribute to the 
propagator (|l6|). This is what one expects, of course, since static fluctuations are excluded from the path integration 
in (§). 

The propagator for the ghosts can be derived in the same way as for the gauge fields. In the Coulomb gauge it is 
equal to the propagator for the temporal fluctuations 

A s c h (i-<',fc) = A° c °(i-i',fc). (22) 

The important conclusion is that the IR-limit of the propagator Ap for the transverse quantum fluctuations is 
0(k~ 1 ) and, therefore less singular than the IR-limit (D(k~ 2 ) of the standard thermal propagator D^. So the severe 
IR divergences that are present in massless thermal field theories in perturbation theory do not occur in a perturbative 
evaluation of the influence action. What we have achieved is a clean separation between quantum fluctuations that 
may be treated perturbatively, and non-perturbative effects contained in the 3D path integral (0). 



D. Thermal weight and equation of motion 

The classical fields are defined on the entire contour C. The classical action S and the effective action W in (|t]), 
contain contributions from classical fields on the entire contour. The contributions from C3 and C\i differ completely: 
the former contributes to the thermal weight whereas the latter determines the time evolution. This difference may 
be made more explicit by demanding 

A cl {t, x) = A cX (t - ia, x) , A° cl (t, x) = -A° cl (t - ia, x) , (23) 

with t € C\ and t — ia E C^. Then (||) implies E c \(t,x) = ~E c \(t — ia,x). We recall that the temporal gauge field 
and the electric field need not be periodic on the entire contour, so the minus sign in (^) does not cause problems. 
The consequence is that the contributions to the classical action from the forward and backward time-branch cancel 
and the contribution from the imaginary time-branch C3 remains 

iS c i = -/3H[A, £] + f3A°K , (24) 

with the abbreviation A°TZ = -2tr / d 3 xA°R[A,£} 7 and 

tSgh[c c i,c c i,yl c i] = iS gh [C,C,A] , (25) 

for the ghost action. The effective action may be split into a part containing only contributions from fields on C3 and 
a part containing contributions from fields both on C12 and C3 

W[A ch E ch cd, c c i, J] = Wbr[A, £,C,C] + W RT [A ch E ch c ch c ch J] (26) 

with Wdr the influence action of static dimensional reduction and the real-time influence action VKrt- The different 
roles plaid by these two actions becomes clear by noting that Wdr is source independent, while a consequence of j2^ ) 
is that the real-time action vanishes at zero source Wrt| j=0 = 0> as f° r sca l ar theories |l8||. Hence, the influence action 
of static dimensional reduction gives the quantum corrections to the thermal weight, while the real-time influence 
action modifies the source term. 

Using the identities (p4]), ( p5| ) and (^) for the generating functional (^), we may write 

Z[j] = J DAV£DCDC5(V ■ A) exp (-j3H[A, £} + PA°U + iS gh [C,C, A] + iWbrJ 

exp (iW RT + ij ■ A cl ) , (27) 

where the first exponent is the effective thermal weight and the second exponent is the effective source term. 

We still have to make a choice for A c \ on the real-time part of the contour. The generating functional is independent 
of this choice and the idea is to minimize the quantum corrections. This may be implemented by demanding 

Sj Wr T [-Acl, E ch cd, cci; J] = . (28) 

3=0 
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To rewrite this requirement as an equation of motion for the classical field, we define the effective action 



T[A ch A] = W RT [A ch E ch c ch c cl ; J]- J- A, 5jW KT = A (29) 

where A is a general fields on C12 (not necessarily equal on C\ and C2). The effective action depends on the zero 
modes A, £,C,C (not explicitly shown) and the field A on C12. Given a source J the field A satisfies 

5 A Y[A ch A]\ A=A = -J . (30) 

Imposing ( p9| ) this gives an equation for the classical field on C12 

S A T[A ch A]\ A=0 = -d A S 12 \ cl . (31) 

We introduce T[A] := T[A,A = 0], where we allow for a general field A on C12. Using the identity ^rfAd, t4]|^_q — 
(5^r[^4]| cl (which follows immediately from a generalization of (6.4) in [[19]), we may write ( |3l| ) as 

MS 12 [A]+r[A])U = 0. (32) 

The action T contains the loop corrections to the classical equations of motion. These loop corrections are calculated 
with the propagator A on the internal lines. 

Note that we may set j = in (|2§| ) and ([32J) , since we have initial conditions on the classical fields at a finite time 



tin and therefore (32) is not a trivial requirement. This in contrast with situations where one has the asymptotic 
condition A c \ — > 0, and it is necessary to keep the source in order to obtain non-trivial solutions. 

A useful property of the vertex functions defined by T[A] is that they satisfy the Slavnov- Taylor identities on C12. 
This can be seen by noting that setting the classical fields equal to zero on C\i gives the standard action for the gauge 
fields ( |l3| ) and the ghost fields. Therefore the action is invariant under BRS transformations of the quantum fields A. 
Also the integration measure DAS (VA) DcS (Vc) DcS (Vc) is invariant under a BRS transformation on C12, since V 
does not act on C12. From the BRS invariance the Slavnov- Taylor identities follow in a standard manner. 



III. QUANTUM CORRECTIONS 
A. Thermal weight 

In this section we use the previously derived Feynman rules to calculate the leading-order quantum corrections to 
the thermal weight in the generating functional (p7|). 

Let us assume we may write the effective action as an infinite sum of local terms, symbolically 

iW D n = 2/3 tr y d 3 xJ2K(d) p (A) Q (£) R (C) s (C) s , (33) 

where the ghost and anti-ghost fields have the same power since these are Grassmann fields. The factor (3 arises from 
the imaginary time integration over C3 . The above equation is a derivative expansion and is valid if the momenta of 
the fields are small k « T . By power counting and a one-loop estimate we find for the coefficient 

K ~ g Q+R+*s T i-N ^ (34) 

with N = P — Q — 2R — 25. We see that terms with 4 — N < are suppressed for small momenta k << T and small 
fields gA, gC « T, g£ << T 2 Note also that the superficial divergence structure is given by the same power 
counting: K is superficially finite if 4 — N < 0. 

We are interested in the leading quantum corrections and we concentrate on terms with 4 — N > 0. Since the theory 
is invariant under static BRS transformations the relevant terms in the effective action must have the form 



iW DR [A,£,C,C] = 2/3tr J d 3 



-^KxT* 1 ? 1 * 1 + ^K 2 £ l £ l - K 3 £ l V l A° + (V l A°) (V l A°) (35) 
-\ml(A ) 2 + ^A (^ ) 4 + K 5 (Cd k d k C - gC[A k , d k C]) + . 



with T kl and T> 1 the static field strength and covariant derivative, respectively. This is the Hamiltonian form of the 
effective action as found in . 
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The coefficients K\...K^^m^ and Ao in ( |36| ) may be calculated in perturbation theory, since the subtractions in 
(|l8|) , (|2(i|) and ( p2|) imply that the perturbative expansion is not plagued with the severe IR divergences of thermal 
field theory. 

The mass for the temporal gauge field equals the well known result for the Debye mass 

ml = -2. 9 2 iV^ J H^nflfc) = l^NT* , (36) 

with k = We use dimensional regularization to render the divergent (spatial) momentum integrals finite; d = 3— 2e, 
fi is a reference mass. The distribution function of the quantum fluctuations is n(k) — n(k) — T/k, and comes from 
the thermal part of the propagator (^p|). The subtraction T/k introduces a linear divergence in the mass (|3^). 
This divergence is temperature dependent in contrast to the 4D divergences of quantum field theory, and acts as a 
counterterm for the same divergence arising in the dimensionally reduced theory. In dimensional regularization power 
divergences are set to zero and the subtraction does not contribute to the mass. Hence the result is the leading order 
contribution to the Debye mass. 

The other constants appearing in (|3^) are dimensionless 4 — N = 0. We find 

* = *> = /!ff< i+ ' Mk)) - = (H- ct • 

A , = WW _^(1 + |)-,,, ( 37 ) 

with | = i + 7 — log 4 ^ and 7 Eulers constant and c.t. the appropriate counterterms. Note that the subtraction 
in h avoids a linear infrared divergence in K2, K3 and K4. 

The constants K\ , Kr and Ao have been calculated before in the Lagrangian formulation of dimensional reduction 
and may be found in [ppCfl for covariant gauges. Integrating out the static momenta £ yields the dimensionally 
reduced action in the Lagrangian formulation. 

Unlike for the linear divergence in d36|), temperature- independent 4D counterterms are sufficient to render the 
constants K\, K$, Ao finite. The constants Ki,...,Kk are subleading in g compared to the classical action. The 
term Xq(A°) 4: is small compared to the term m§(A ) 2 for not too large fields A° ~ T, since Ao ~ g 4 ; see (|34|). Hence, 
we conclude that only the m^A ) 2 is not perturbatively small compared to the classical action. This is the only 
static HTL and has to be included in the thermal weight for the effective classical theory, while the other subleading 
contributions can be treated as perturbations. 



B. Hard thermal loop contributions to the equation of motion 



We now turn to the equation of motion ( p2| ) . The effective action T contains all quantum effects and is a complicated 
expression. However, since we are interested in low- momentum correlation functions k << T, we may consider only 
the leading-order terms to the equation of motion at low momenta and weak-coupling. 

The leading contributions to Y at low momenta are the well-known HTL's. We review the order estimates [[T|]2~l|l 
for the HTL contributions to the action T and discuss the modifications due to the subtraction in the quantum 
propagator. 

In the effective classical theory one has two length scales (gT)~ l and (g 2 T) -1 . These two scales correspond to 
two different excitations; the plasmons that are screened over a length scale (gT)^ 1 [ p2| , and the attenuating modes 
p3[ that are screened non-perturbatively over a length scale (<7 2 T) _1 due to a magnetic screening mass in 3D gauge 
theories 0. 

Let us first consider the momentum and energy scale k ~ gT, fco ~ gT of the plasmon modes. We assume that 
the classical fields are of the same order as the initial fields A c \ ~ A. An estimate for the initial fields is found by 
requiring that the fields are not suppressed by the thermal weight: (3H ~ 1 — > A ~ ^fgl '. This yields the following 
estimates for terms in the HTL action 

r^S~ ff 2 +^T 4 , (38) 

where L^ TL are the retarded vertex functions, as they appear in the equations of motion. By definition of the HTL's, 
the classical and HTL contributions are of the same order. The interaction terms are small compared to the quadratic 



parts and therefore (resummed) perturbation theory is valid for the plasmon scale 24 2 
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Secondly we consider the momentum scale k ~ g 2 T and energy scale fco ~ <? 4 T of the attenuating modes EJ. 
These are only present in the spatial gauge fields, since the temporal gauge field is electrically screened with Debye 
screening length (gT)~ 1 . Again we assume that the classical fields and initial fields are of the same order, yielding 
A c \ ~ A ~ gT. We have the following estimates for terms in the classical and HTL action 

rP TL A% ~ g 2 - n T"- n k Q A^ ~ g 6 T 4 . (39) 

Again the HTL contributions are as large as the classical contributions. We see that higher order terms are not small 
compared to the quadratic part of the action, which means perturbation theory fails. It will not be possible to expand 
in the non-linearities of the effective classical theory. The estimate for the classical and HTL action ([39]) is of the 
order where the perturbation expansion breaks down for the free energy jj) . 

In thermal field theories large contributions to vertex functions may also come from the IR-region H. These may 
be as large as the HTL and classical contributions. However these contributions are not present in (|32|) since we 
use the subtracted propagator, which implies that IR contributions are suppressed by a factor /3k ~ g. The large 
IR contributions are included in the effective classical theory. This leads to the conclusion that the HTL's have to 
be included in the effective equations of motion and that at the scale g 2 T a perturbation expansion for the effective 
classical theory fails |ll],[l3| . 

Let us consider HTL's of the quantum fluctuations. These HTL's are calculated with the Feynman rules given 



in section QC. We focus on the retarded vertex functions, since these appear in the equations of motion. We may 
restrict ourselves to vertex functions with all (external) times on Cn , since also correlations between the initial fields 
and the classical fields at later times can be expressed in these vertex functions as explained in [p78[] . Since the HTL's 
consist of one loop diagrams we only need the interactions and quantum propagator on C12. The interactions are 
the usual 3- and 4-point interactions, as we see from (|l3|). The only difference with the usual Feynman rules of 
thermal field theory is in the propagator, where the Bose-Einstein distribution function is replaced by a subtracted 
distribution function in the quantum propagator ubeW) — > n(k) — nfls(fc) — T/k. However, we shall now argue that 
this does not affect the well-know HTL expressions |2l]] in dimensional regularization. Retarded one-loop diagrams 
contain one distribution function plfl . The momentum integrals of the usual HTL's have a cut-off at k ~ T due to 
the Bosc-Einstcin distribution function yielding a T 2 behavior. The subtraction in the distribution function of the 
quantum fluctuations introduces a linear divergence in the HTL's. Since the classical linear divergences arise only 
at one- loop and the generating functional is finite (including 4D counterterms) the linear divergence in the HTL's 
acts as a (non-local) counterterm for the classical theory. Since we use dimensional regularization, which sets linear 
divergences equal to zero, the HTL's of the quantum fluctuations are just the usual HTL's. In appendix A, we analyze 
the linear divergences using a momentum cut-off on the momenta of the initial fields. It turns out that the linear 
divergences enter in the plasmon frequency only. Such a linear cut-off dependence has been proposed before from a 
matching argument |l3|,^6| . The analysis here shows how such a counterterm arises from integrating out the quantum 
fluctuations. 

The HTL's we thus obtain are gauge invariant, just as the usual HTL's. This is as a consequence of BRS invariance 
of r and the HTL approximations [Q . 

IV. EFFECTIVE CLASSICAL STATISTICAL THEORY 
A. Exact generating functional 



In section [I D an effective classical statistical theory was defined. Following the literature |f|, |f Qj , P-3| , it was argued 
that the HTL contributions should be included in the effective Hamiltonian and the equation of motion. We argued 
also that the HTL's of the quantum fluctuations equal the usual HTL's in dimensional regularization. With this the 
exact generating functional can be rewritten as 

Z[j] = J DABSDCDCS (v • A) exp-[3H cS {A,£] + iS gh [C,C,A] + iWgg 

expiWKr + ij-Ad , (40) 

where we have defined the effective Hamiltonian 

H cS = H- PA°K + 2tr J d 3 x^m 2 (A ) 2 , (41) 
and the static non-HTL corrections 
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WSS = W m . + 2z/3tr J d 3 x^ml(A ) 2 . (42) 

In appendix B, it is shown this Hamiltonian corresponds to the conserved energy of the classical subsystem. 
The HTL equation of motion is 

— v,J^ dlfU*(x,J)v u F%'(x')e-« t -V , (43) 

with = g 2 NT 2 /9 the plasmon frequency. The angular integration is over the direction of u, |u| = 1. Furthermore, 
x 1 = (t',x — v(t— t')) and the parallel transporter U c \{x,x') — Pexp(— ig J^dz^A^z)) , with 7 a straight line from x 
to x' . The damping exponential follows from the e-prescription of thermal field theory [l7| . 

The HTL equation of motion has been derived by Blaizot and Iancu p7[ , using a kinetic approach. They have 
shown the HTL equation of motion is a linearized Vlasov equation (linear in the deviation from equilibrium), with 



the r.h.s. of (43) an induced source. 



The initial conditions for the equation of motion fl43| ) are 

A k cl (t in , x ) - A k (x) ; F° k (t in ,x!)=£ k (x). (44) 

at a finite time t m . The Lagrange multiplier A is not a dynamical field, so A does not act as an initial condition. 
The path integration in (f40|) accounts for the thermal average over the initial conditions (E3). 



B. Classical approximation 

Note that the generating functional (|4C|) is still exact. The necessary approximation to obtain an effective classical 
theory is to neglect W^r and Wrt in p0|), which yields 

Z c i [j] = J DADSDCDCS (v ■ Aj exp -{3H eS [A, S] + iS gh [C, C, A] + ij ■ A ci . (45) 

Let us argue why and Wrt in (0) may be neglected for the calculation of correlation functions at low momenta. 
The action contains only the non-HTL contributions because of (|d]) and ([l2]). Since these are small compared 

to the classical and HTL contributions in H c g the action Wgj^ may be neglected. 

The other approximation Wrt = can be justified as follows. We consider as an example the symmetric transverse 
two-point function 

D T (h-t 2 ,k) = ^(^.-^(^(ta.fc) + (ti «-» t 2 )) , (46) 
which can be obtained from the generating functional 

D T (t 1 - t 2 , k) = -^T lJ [Sj^ti-ia^jjitz) + S ji(ti)^j } (t 2 -ia)] , (47) 

where the a indicates that the source derivative is taken on the backward contour C2, see fig. 1. Since the classical 
fields satisfy the effective classical equation of motion the contributions involving SjWnr may be neglected. We then 
find that the two-point function can be separated in a classical and a quantum contribution 



D T (ti - t 2) k) = D T \t! - ti, k) + D T (h - t 2 , k) (48) 



with 



b^{h-t 2l k) = i<4 1 (*i ! -^W4l(*2,fc))cl , (49) 

D T (ti -t 2 ,k) = -7([<Jj(ti-to-)^j(t 2 ) + fij(t-L)8j(t 2 -i<r)\ W RT )c\ ■ (50) 

The brackets (...) c i denote the classical average over the initial conditions defined by ([f5f). In temporal Fourier space 
these contributions take the form 
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with the classical spectral density 



D$(ko,k) = p c \k ,k)n cl (k ) 



D T (k ,k) 



p c \k a ,k) = 21m 



(fl(ko,k)-[l + 2n(ko)] , 



ik + ie) 2 + k 2 +H c R (k ,k) 



(51) 
(52) 



(53) 



and Hft(ko,k) the classical retarded self-energy. A similar expression can be given for the spectral density p q with 
self-energy 11^. The appearance of n = tire — T/ko i n (|5^ ), instead of the Bose-Einstein distribution function, is a 
consequence of the subtraction in the free propagator fl20|). In fact, in the free case (H^. — 0) the two-point function 
([H|) corresponds to the free propagator fl2C| ) . Note that in the case that the classical fields satisfy the effective equation 
of motion (32) exactly SjWur vanishes identically and as a consequence the self-energies irjj,lT^ are equal to the full 
self-energy. We have required that the classical fields satisfy the HTL equation of motion, in this case the leading 
contribution to both the self-energies 11^,11^ is the HTL self-energy. Indeed, the action Wrt contains the HTL 
self-energy and all other HTL vertices. However the quantum contribution D T to the symmetric transverse two-point 
function can be neglected at small momenta \ko\ <<T, since 



[1 + 2n(fe )] 



0<k Q «T 



(54) 



Therefore we may neglect the action Wrt for the calculation of (symmetric) correlation functions at low momenta. 
Note that the subtraction in the two-point function for the quantum fluctuations is essential to justify the classical 
approximation (^) for the calculation of low-momenta correlation functions. 



V. CONCLUSION 



We have derived an effective classical theory by integrating out quantum fluctuations. The effective equation of 
motion (|32|) is is BRS invariant, which in the HTL approximation implies gauge invariance. The Feynman rules for 
the calculation of the effective equation of motion are the same as usual except for a subtraction in the propagator. At 
real-times (on the contour C12) this subtraction enters only in the distribution function h = nsE—T/k. We would like 
to stress that this subtraction is a direct consequence of the extraction of the path integration over initial fields from 
the quantum path integral. This subtraction is essential in two respects. It allows one to argue that the corrections 
from Wrt to symmetric correlation functions are small and therefore the classical approximation is justified. Secondly, 
the subtraction introduces linear divergences in the HTL equations of motion and the effective Hamiltonian that act 
as counterterms for the linear divergences in the classical theory. We used dimensional regularization which sets linear 
divergences equal to zero. Since a non-perturbative treatment of the classical theory requires a lattice implementation, 
the generalization to a lattice regularization is important, see p8[ . 

There are still logarithmic divergences present in the effective classical theory (|4^). Therefore the classical theory 
( |45"| ) has to be formulated with some cut-off A, which has to be large compared to the plasmon scale gT << A. Since 
logarithmic divergences are suppressed by a factor k/T ~ g,g 2 (at one-loop) or g 2 (at two-loops) the cut-off may be 
taken as large as the temperature A ~ T, since then the requirement (k/T) log(A/<?T), g 2 \og(A/gT) << 1 is fulfilled. 

In principle, T and Wbr furnish the logarithmic counterterms. If these can be extracted explicitly the result is 
a cut-off independent classical approximation. The counterterms for the one-loop divergences of the classical theory 
may be extracted from the one-loop logarithmic divergent contributions in T and Wdr- However, it is unclear how one 
obtains the counterterms for the two-loop logarithmic divergences. Yet, we know that the full one-loop contributions 
together with the two- loop logarithmic divergent contributions should be sufficient. This suggests that an analysis as 
in p3] for static dimensional reduction may be generalized successfully. 
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APPENDIX A: HARD THERMAL LOOPS WITH A SHARP CUT-OFF 



To explicitly see the linear divergence in the HTL's, we use a sharp cut-off A on the momenta of the classical 
fields at t{ n to regularize divergent momentum integrals. This imposes the restriction k < A on the initial fields and 



consequently on the subtraction (|21j) on C\i. Repeating the arguments of section [II B, we see that this reduces to a 
replacement of the Bose-Einstein distribution function of the usual HTL's by h A 

n(k) -> n A (k) = n(k) - -6(A - k) . (Al) 
k 

In the HTL expressions for the retarded vertex functions, we use the expansion 

dkk 2 h A (k +p cos 6) = J dk [k 2 h\(k) - 2k h A (k)pcos 9 + ...] (A2) 

in the external momenta p. The angular integration decouples from the k integration, as usual in the HTL approx- 
imation. Performing the k integration we find the familiar HTL contributions to the vertex functions, but with a 
A-dependent plasmon frequency 

(^) 2 =^i-^9 2 NAT. (A3) 

This le ads to the conclusion that a linear divergence appears in the plasmon frequency. By the same reasoning as in 
section 



III B , it is clear that this linear divergence acts as a counterterm for the classical theory. 



In the literature |T(]|3(],[nJ , often a cut-off is introduced according to 

n(k) -► n A (k) = n(fc)G(A - k). (A4) 

In the classical approximation, one chooses gT « A << T and approximates the soft field by a classical field 
fl0| , [30f . The classical theory then contains the leading order contributions to low- momentum correlation functions up 
to corrections of order /3A. However, the result including all contributions from integrating out the hard modes is not 
exact anymore. The reason is that an approximation is made on the low-momentum field, namely that it behaves 
classically. This implies that the dependence on the cut-off docs not cancel out in the final result. A linear dependence 
on the cut-off will remain. 



APPENDIX B: EFFECTIVE HAMILTONIAN 



In this appendix we consider the classical subsystem, that is the generating (E^) functional without the quantum 
corrections — = Wrt- We show that the effective Hamiltonian in (|41|) corresponds to the energy of the 

classical subsystem. We follow Blazoit and Iancu [^7| and introduce an auxiliary field W°(x, v), that is also a function 
of the direction of v; v — (l,v) and \v\ = 1. This allows for a local formulation of the equations of motion ( [43|) 



K,<](z)=3u, 2 



pl 



47T 



v„W°(x,v) 



[v„D» h W°(x,v)] =F$(x)v p 



with the covariant derivative D'f, = <9 M + igA't,. The conserved energy of this system (B1,B2) is 27 



E = -tr / d'r 
In particular, at t- ln the energy is 



dfi. 



E 2 {x) +B 2 (x) +3^ I — W°{x,v)W°{x,v) 



d 3 x +B 2 {x) + 3^ / W°{t in ,x,v)W°{t in ,x,v) 



(Bl) 
(B2) 

(B3) 
(B4) 



The value of the auxilary field at t- lTL follows from the prescription to keep the fields constant in time before t m in 
(p3|). In a time-independent background field the vector current is zero , hence, the vector current at t m vanishes. It 
follows that W°(t in , x, v) is independent of the direction of v. A calculation yields 
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W°(t in ,x,v) = A°(t in ,x) . 
Using the equation of motion at ti n , we have 

[D l ch F$] (t ia , x) = U{x) = 3w^°(tto, x) 



Inserting (B5) and (B6) in (B4), we find for the energy of the classical subsystem 



E : 



-tr 



d 3 a 



P(x) +B 2 {x) + —^-TZ 2 {x) 



Hi 



(B5) 



(B6) 



(B7) 



Let us consider the effective Hamiltonian Since we neglect the influence actions in (|l|), we can simply integrate 
out the zero mode of the temporal gauge field A . This yields the effective Hamiltonian H c h = E in the thermal 
weight. 

An effective classical theory based on the equations of motion (B1,B2), has been proposed by Iancu The 
average over the in itia l conditions is taken by a path integration over A,£ and W°(t in , x, v ), constrained by Gauss' 
law. We see from (B5), that restricting the path integration over ^-independent auxiliary fields at t in , is consistent 



with the inclusion of the correlations between initial fields and the fields at lat er t imes according to the prescription 
in (^3|). The value for initial auxiliary field is then determined by Gauss' law (|B6| ). This yields a somewhat simpler 
effective classical theory than proposed by Iancu, since no path integration over W (t- ln , x, v) has to be performed. 
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FIG. 1. The time contour C. 
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